Spectrum and diffusion for a class of tight-binding models on hypercubes 



OO 
On 
On 

> 
o 

m 

(N 



-a 
c 

o 
o 



> 
m 



oo 

ON 

c3 



I 

c 

o 
o 



X 



Julien Vidal 1 , Remy Mosseri 1 and Jean Bellissard 2,3 
1 Groupe de Physique des Solides, Universites P. et M. Curie P6 et D. Diderot P7, 
Tour 23, 2 place Jussieu, 75251 Paris Cedex 05 
2 Institut Universitaire de France 
! Laboratoire de Physique Quantique, UMR5626 associee au CNRS, IRSAMC, Universite Paul Sabatier, 
118, route de Narbonne, F-31062 Toulouse Cedex 4, France 

We propose a class of exactly solvable anisotropic tight-binding models on an infinite-dimensional 
hypercube. The energy spectrum is analytically computed and is shown to be fractal and/or ab- 
solutely continuous according to the value hopping parameters. In both cases, the spectral and 
diffusion exponents are derived. The main result is that, even if the spectrum is absolutely contin- 
uous, the diffusion exponent for the wave packet may be anything between and 1 depending upon 
the class of models. 
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The interplay between the energy spectrum and the 
diffusion process for quantum systems is still an open 
question. Indeed, recent studies have shown that the 
spreading of a wavepacket is determined by several ex- 
ponents that only depend on the nature of the spectrum 
(absolutely continuous, singular continuous, pure point, 
or any mixture) [|l]-^). In particular, for quasiperiodic 
systems, such as quasicrystals (QC), it has been shown 
that an anomalous behaviour of the wavepacket spread- 
ing occurs when the spectrum is singular continuous Q 
(see below). 

In such systems, scaling laws occur at many levels. We 
will concentrate upon two classes of scaling exponents 
due to their importance for electronic and transport 
properties. The first one concerns spectral exponents, de- 
fined as jjH+gz dN{E') ~ 5E a W as SE -> 0, where M 
is the local density of states (LDOS). An absolutely con- 
tinuous spectrum in some interval / of energy implies 
a{E) = 1 for E E I and E in the spectrum. A pure 
point spectrum in /, namely the LDOS is a sum of Dirac 
peaks in /, implies a(E) = on the part of the spec- 
trum contained in I. Finally, if < a(E) < 1 for E 
in some part of the spectrum implies singular continu- 
ous spectrum there. This is actually what happens for 
ID QC HD such as the Fibonacci chain. Such singu- 
lar continuous spectrum occurs also for ID chains with 
a potential given by a substitution sequence [fz| — |l0|] ■ The 
problem of 2D electrons on a square lattice under mag- 
netic field that can be mapped onto a ID Harper equation 
JTT| also display a singular continuous spectrum jl2j for 
incommensurate fluxes, and the spectral exponents have 
been computed numerically p3| . A number of non rigor- 
ous results suggesting singular spectra were obtained for 
models on higher dimensional QC such as the labyrinth 
model PI , the octagonal tiling [flEf , in a regime of small 
hopping. At last let us also mention recent studies on 
Jacobi matrices associated to iterated functions systems 
(IFS) QUI generalizing Julia sets pf^ leading to the 



calculation of the spectral exponents. 
The next class concerns diffusion exponents, given by 
LE{t) ~ t@( E ^ as t — * oo where Ls{t) characterizes 
the spreading of a typical wave packet of energy E af- 
ter time t J2l| (see eq. ([l6]) for a definition of the expo- 
nents). Ballistic motion such as the electronic motion 
in a perfect crystal, corresponds to f3 = 1. Strong local- 
ization is defined by demanding that sup t>0 L^it) < oo 
implying (3(E) = 0. The weak localization regime appear- 
ing in weakly disordered metals generally corresponds to 
j3 = 1/2, namely to the case for which the quantum evo- 
lution mimicks a classical diffusion such as a Brownian 
motion. Such a behaviour of the quantum diffusion also 
occurs in the Random Phase Approximation (RPA) [£lj . 
If (3 takes on values different from 0, 1/2, 1 we will speak 
of anomalous diffusion. In QC's for instance, numerical 
results show that (3 may be anything between and 1 
and is model dependent For realistic QC's, such 

as i-AlCuCo, an ah initio calculation leads to (3 — .375 
at the Fermi level Q. 

Spectral and diffusion exponents are related through 
the Guarneri inequality B, 



(3(E) > 



<E) 



(1) 



where d is the dimension of the system. This inequality 
implies that for absolutely continuous spectrum, ballis- 
tic motion always occurs in ID (see however ]23 , B4]| ), 
whereas in higher dimension (3 > 1 / d. One of the main 
question addressed recently in this respect is whether f3 
can be computed more directly from a. It seems that this 
is indeed the case for Jacobi matrices |17||l8fi namely ID 
chains with nearest neighbour interactions. One impor- 
tant result of this paper is precisely to show the opposite 
in the extreme case for which d = oo, namely that there 
is no relation whatsoever between the a's and the /3's. 
This gives a negative answer to the question raised by 
Lebowitz (see 0). This is because the spectral expo- 
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nents characterize only the spectral measure (the LDOS) 
of the Hamiltonian H, independently of any other type of 
observables. On the other hand, the diffusion exponents 
involves the interplay between the Hamiltonian (through 
the quantum evolution) and the position operator X, or 
even better, the current J — i[H 1 X]/h. The link between 
the diffusion exponent (3 and the pair (H, J) is still not 
precisely established, even though it has been related to 
the current- current correlation function |p5[pl| . How- 
ever, Jacobi matrices are very special since the position 
operator is defined by mean of the orthogonal polyno- 
mials associated to the spectral measure (the LDOS), so 
that (3 is defined through purely spectral properties. 

In this letter, we consider a family of anisotropic tight- 
binding models in an infinite dimensional hypercubic 
structure and show that depending upon the explicit 
form of the hopping parameters, it is possible to shift 
from an absolutely continuous spectrum to a singular 
continuous spectrum. Moreover, we are able to adjust 
the fractal dimension of this spectrum, fine-tuning a sin- 
gle parameter that drives the transition. In addition, we 
show that depending on the hopping term law, one can 
face an absolutely continuous spectrum and a somewhat 
anomalous diffusion for which the mean square deplace- 
ment L(t) can either scale as log t or as t a with < a < 1. 
We first introduce some mathematical tools that are use- 
ful for a careful analysis of the structure we are dealing 
with. We then characterize the energy spectrum for dif- 
ferent types of tight-binding hamiltonians and discuss the 
nature of their spectral measures. Finally, we compute 
the autocorrelation function and the mean square dis- 
placement of a wavepacket for the different models. 

A (i-dimcnsional hypercube A^ is the set of vertices of 
a cube of size I in an e?-dimensional space. The infinite 
dimensional hypercube A is defined by: A = lj d>1 A^. It 
can therefore be seen as the set of sequences e = (£k)kLo 
where £ {0, 1} and Sk — for all but a finite number 
of fc's. We endow the set {0, 1} with the group structure 
given by the addition modulo 2, so that A becomes a dis- 
crete countable group for the coordinatewise addition. It 
is also convenient to introduce its dual group B (B stands 
for Brillouin) , which is the counterpart of the quasimo- 
mentum space in a perfect crystal. B can be described as 
the set of all sequences a — (crfc)fco wrtn a k = ±1- B is a 
compact abelian group with the pointwise multiplication 
and the product topology. The duality between A and B 
is given by the characters: 

oo 

V<7eB,VeeA, x*(£) = H (T k- ^ 

k=0 

In this formula, the product is finite by construction. 
Moreover, x<y{e + e') = Xo{z)x<t{z') and Xv+a'i?) = 
X(j( £ )x<j'( £ ) ■ These characters play the role of the Bloch 
phase exp (ika) in a crystal, where a is a period of the 



translation group, and fc is a quasimomentum. While on 
A the Haar measure is the counting one, the integral of 
a continuous function / on B is defined as: 

f d*f(*)=Xm*r E /<*)• (3) 

The Hilbert space of physical states is TL = £ 2 (A), 
namely the set of sequences -0(e) indexed by A (the wave 
functions), such that: 

ihi 2 = Ei^ £ )i 2 <+°°- ( 4 ) 

eSA 

A canonical orthonormal basis is provided by the states 
|e > vanishing everywhere but on the "site" e. The 
Fourier transform of a wave function ip € H is the func- 
tion on B formally defined by: 

%) = E^W s )' ( 5 ) 

EfEA 

This function actually belongs to L 2 (B), namely it is 
square integrable on B (with respect to the Haar mea- 
sure) and the Parseval identity holds true namely | l^l | 2 = 
l-^^ll 2 = Jb da\!Fip(a)\ 2 . Therefore we get two unitarily 
equivalent representations of the Hilbert space of states. 

The translation operators T(a), (a £ A) are acting on 
TC as follows: 

T(o)^(e) = tl>(e - a) . (6) 

Equivalently, T(a)\e >— \e — a >. Note that a — —a in A 
due to the addition modulo 2, so that T(a) 2 — 1, Va 6 A. 
In addition T(a) — T(a)^ as can be easily checked, so 
that there is a infinite set of mutually commuting uni- 
tary and self-adjoint operators. The spectrum of such 
operators is made of two eigenvalues (with infinite multi- 
plicities) namely ±1. Through a Fourier transform, T(a) 
becomes the operator of multiplication by Xa{o)- Parti- 
cularly, if a — efc, where is the sequence in A with all 
coordinates vanishing except the fc-th one, = T(efc) 
becomes simply the operator of multiplication by <jfc. 

We consider the following class of tight-binding hamil- 
tonians on Ti: 

oo 

H = E tkTk ■ ( ? ) 

In order that H be self-adjoint we need t^ 6 R. By a 
simple unitary tranformation, one can choose tk > 0. 
The coefficient tk denote the "transfer" or "hopping" 
term in the fc-th direction. H is bounded if and only 
if X^fc < +oo- It is self-adjoint (but not necessarily 
bounded) if ^2 1 2 < +oo. In what follows, we will assume 
that this latter condition holds. By Fourier tranform, H 
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becomes the operator of multiplication by E(a) where E 
is called the band function and is given by: 



depend upon the explicit choice of ip as it turns out. Us- 
ing the previous relations, one gets for any E: 



e(o-) = y^o-fctfc 



(8) 



fe=0 



This function is real and square integrable on B with L 
norm: 



daE{af = 1 



(9) 



k=0 



The spectrum of H (resp. its spectral measure), is then 
given by the image of B in R (resp. of the measure da) 
under the function E. Note that if H is bounded E is 
continuous with = sup aeB \E(a)\ = J2kLo t k- 

Spectral properties can be studied through the auto- 
correlation function: 



P(s) 



< 0\e lsH \0 > 



dn{E)e l 



R 



(10) 



where |0 > denotes an origin site where we initially local- 
ize a wavepacket and /j, the corresponding spectral mea- 
sure. Note that the translation invariance of H allows 
us to choose any site of A as initial condition. If P is 
integrable over R, then /x is absolutely continuous (the 
converse may not be true). One can alternatively use the 
temporal correlation function: 



i r* 

C(t) = - J dsP( 



s) 



(11) 



that is the time-averaged version of P. The spectral mea- 
sure is purely continuous (singular or absolutely continu- 
ous) if and only if C(t) — > as t — > oo (Wiener criterion) . 
An elementary computation using (|^) leads to: 



P(s) = f[ cos 2 (st k ) . 



(12) 



k=Q 



This infinite product converges since ^t 2 . < +oo. 
We define the position operator as follows. For k G N, 
Xk denotes the operator of multiplication by Sk in Ti. It 
commutes with 7} for I ^ k whereas T^X^T^ 1 = 1 — Xk, 
and since Tj* = 1 it follows that: 



cos s + iTk sin s 



(13) 



Let us define the mean square displacement by: 

oc 

l 2 e ( S ) = y,<^( x ^- x ^ 2 ^> > ( 14 ) 



fc=0 



where Xk(s) — exp(isH)Xkex-p(—isH) denotes the 
Heisenberg representation of Xk and \ip > is an initial 
state with energy close to E. This expression does not 



L 2 E (s)=L 2 (s) 



] sin 2 (st k ) 



(15) 



fe=0 



The diffusion exponent (3 is given by L(s) ~ as s — > oo. 
It does not depend on E. A rigourous way to define a 
power law asymptotic behaviour is given as follows (see 
|^l| , p6| for more details): a function / of a real variable 
s behaves as when s — * oo if: 



ds 

c l+6 



f(s) 



c> 



(16) 



converges for b > [3 and diverges for b < (3. In addition, 
if the function / can be written as a series: 



f(s) = J2F(st k ) 



(17) 



fe=0 



where F is a positive bounded real function, such that 
F(x) = 0(x 2 ) for x ~ 0, and (tfc)feeN a set of positive 
number such that X)fc°=o ^fc ^ °°j then the exponent (3 is 
given by: 



(3 = inf {b € R 4 



fc=0 



4 < oo} 



(18) 



It is clear that this definition is particularly convenient 
for our purpose since L 2 is exactly of the form (^) . 
The closed forms obtained for the three observables previ- 
ously discusssed (energy, autocorrelation function, mean 
square displacement), allows us to study the spectrum 
and the quantum diffusion for various classes of hopping 
terms. 

A first interesting class of models consists in choosing 
an algebraic scaling of the hopping parameters tk ~ k 1 , 
namely lim^oo fc 7 ifc = t, with 7 > 1/2. In this case, the 
spectrum is bounded if 7 > 1 whereas it is unbounded if 
7 < 1. Moreover, one can prove (see Appendix) that 



P(s) < cie" 



,1/7 



(19) 



where ci,C2 are two positive constants. This shows that 
the spectral measure is always absolutely continuous and 
also infinitely differentiable. This also implies that the 
correlation function decays as 1/i. In addition, according 
to expression (|l8|), it is obvious that: 



L 2 (s) 



a/7 



(20) 



Hence, the diffusion exponent is j3 = I/27 which can take 
any value in ]0, 1[ even though the spectrum is always ab- 
solutely continuous. 

Another interesting case is tk = (q — l)/V c+1 - ) (geomet- 
rical scaling) with q > 1, for which ||i?|| = 1. 
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(i) For 1 < q < 2, the spectrum is nothing but the q— adic 
decomposition of real numbers in the interval [— l.+l]. 
It is therefore gapless and absolutely continuous. 

(ii) For q > 2, the image of E is a Cantor set of zero 
Lebesgue measure, constructed by removing the central 
interval of width 2(1 — 2/q) in the interval [—1, +1] and 
repeating the operation on each of the intervals left. The 
spectrum is a monofractal set with a Hausdorf dimension 
Dh = In 2/ In g pTj]. The spectral measure is the Cantor 
one and gives the same weight to each subinterval. Note 
that the classical tryadic Cantor set is obtained for q = 3. 
For such Cantor spectra, it is shown in ref. [Q that the 
temporal correlation function decays as C(t) ~ t~ D ' 2 , 
where D 2 is the correlation dimension of the spectral 
measure (i. e. of the local density of states). In this ex- 
ample, one has: D 2 = 1 for 1 < q < 2, since the spectrum 
is absolutely continuous, and D 2 = Dh = log 2/ log q for 
q > 2 since the spectrum is then a monofractal set. It is 
important to consider C because the behaviour of P is 
much more complex. In particular, P is sensitive to the 
nature of q. Indeed, it is shown in p8| that: 



Appendix: proof of eq. (|19|) 



Let us consider K > large enough so that t/2k 7 < 
tk < It/k 1 for k > K. Then choose so > large enough 
so that s a t/K~< < vr/2 < s t/(K - l) 7 . For s > s let 
K 1 > K be such that st/2Kj < tt/2 < st/2(K 1 - l) 7 . 
Then 



lnP(s) < ^2 lncos 2 st/2fc 7 , 

k=Ki 

If one sets Xk = k (2/st) 1 ^ 1 , the right hand side is domi- 
nated by an integral of the form 



lnP(s) < (ts/2) 1 ^ / da;lncos 2 (l/x 7 ) 



for s > sq, leading to eq. (UM 



lim P(s) = & q i S \ {2} , 



(21) 



where S denotes the set of algebraic integer numbers de- 
fined by Pisot and Vijayaraghavan p9 30 1. 
Note that the functional relation: 

Vg, P(qs) = cos 2 {s(q-l))P(s) , (22) 

allows us to exactly determine P for q = 2 since: 

P(2s) = cos 2 (s)P(s) ^ P(s) = sin 2 (s)/s 2 . (23) 

According to the identity (|l8|), one obtains a diffusion 
exponent (3 = for any q > 1 whereas the spectrum 
can be either absolutely continous or singular continuous. 
In addition, one can show that L 2 (s) ~ ln(s), with a 
criterion similar to the one given in eq. (fl6|). 

To conclude, these toy models defined on infinite di- 
mensional hypercubes allows us to carefully analyze the 
possible relationships between the spectral measure and 
the diffusion exponents. The first class of hamiltonians 
(algebraic scaling of the hopping terms), shows that it is 
possible to face an absolutely continuous spectrum and 
an anomalous diffusion with a f3 exponent that can take 
any value between and 1. On the other hand, the sec- 
ond case (geometrical scaling of the hopping parameters), 
displays a zero (3 exponent whereas the spectrum can 
be either absolutely or singular continuous. Finally, we 
emphasize upon the importance of these exponents in 
transport properties, especially in quasicrystals, where 
they should be responsible for the anomalous behaviour 
of their conductivity. 
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